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Abstract. We construct the solutions of slowly rotating gravastars with a thin shell. 
In the zero-rotation limit, we consider the gravastar composed of a de Sitter core, a 
thin shell, and Schwarzschild exterior spacetime. The rotational effects are treated as 
small axisymmetric and stationary perturbations. The perturbed internal and external 
spacetimes are matched with a uniformly rotating thin shell. We assume that the 
angular velocity of the thin shell, fl, is much smaller than the Keplerian frequency 
of the nonrotating gravastar, flfc. The solutions within an accuracy up to the second 
order of are obtained. The thin shell matter is assumed to be described by a 

perfect fluid and to satisfy the dominant energy condition in the zero-rotation limit. In 
this study, we assume that the equation of state for perturbations is the same as that of 
the unperturbed solution. The spherically symmetric component of the energy density 
perturbations, Jctqj is assumed to vanish independently of the rotation rate. Based on 
these assumptions, we obtain many numerical solutions and investigate properties of 
the rotational corrections to the structure of the thin shell gravastar. 
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It is now widely believed that the black hole is a possible hnal state of massive stars 
after the gravitational collapse. Several tens of stellar-mass black hole candidates have 
been identihed in our Galaxy. With a rough estimate, then, 10® of stellar-mass black 
holes are expected to exist in our Galaxy. A supermassive black hole of 10® to 10® solar 
masses is also believed to be in the centers of most galaxies. Thus, it is likely that there 
is a huge number of black holes in the Universe. If quantum effects on gravity are taken 
into account, however, it could be that the gravitational collapse of objects comes to a 
halt and furthermore no event horizon forms [1]. If this is the case, there is no black 
hole in the Universe. Some people have seriously considered this possibility and have 
proposed alternative models to the black hole that are almost as compact as black holes 
but have no event horizon. 

The first serious model of such black hole-like objects is the so-called gravastars, 
gravitational vacuum stars, proposed by Mazur and Mottola [2]. During the formation 
of the gravastar, just before the event horizon formation, a quantum vacuum phase 
transition is assumed to occur around the location where the event horizon would have 
been expected to form. As a result, the black hole event horizon formation is prevented. 
Actual equilibrium models of the gravastar are composed of a de Sitter core (a positive 
cosmological constant matter core) surrounded by a relatively thin but hnite-thickness 
shell whose equation of state is given by p = a, where a and p are the total energy 
density and pressure of the shell matter, respectively. By gluing these constituent parts, 
they consider objects having no event horizon nor spacetime singularity. Mazur and 
Mottola also showed that their gravastars are thermodynamically stable unlike in the 
case of the black hole. 

A simplihed version of the gravastar model has been considered by Visser and 
Wiltshire [3]. This simplihed model is called the thin shell gravastar because the 
relatively thin but hnite-thickness shell of the original model is replaced by the single 
inhnitesimally thin shell. They have further made detailed investigations of the 
dynamical properties of the thin shell gravastar and shown that there are thin shell 
gravastars stable against spherically symmetric perturbations. Some other models of 
the gravstar have been considered by many authors, see, e.g., HI El El [7]. 

In standard astronomical observations, the black hole candidates are basically 
chosen either because their estimated mass is greater than some theoretical limit, or 
on the basis of their strong external gravitational held. According to these astronomical 
criteria of the black hole candidate, therefore, all the black hole candidates could be 
alternative models whose mass and compactness are similar to those of the genuine black 
hole if the existence of such alternative models to the black hole is allowed theoretically. 
However, present and planned gravitational wave observatories will give us data with 
which one may test whether black hole candidates are genuine or not and check what is 
likely to be as an alternative model if the black hole candidates are not genuine black 
holes. 
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One of such data is the characteristic oscillation frequency of the black hole-like 
object, which is called the quasinormal mode frequency because it is given by a complex 
number, whose imaginary part represents damping rate of the wave due to the energy 
loss. For the case of gravastars with a thick shell, Chirenti and Rezzolla [6] have 
calculated quasinormal modes of axial parity perturbations. They showed that the 
decay rate of the gravastar quasinormal mode is substantially different from that of the 
black hole when the gravastar and back hole have the same gravitational mass. They 
also found no unstable eigenmode in the gravastars they considered. For the case of 
the thin shell gravastar, Pani et al. [8] have calculated the quasinormal modes of the 
axial and polar parity perturbations. Their main results are the following: (i) The 
quasinormal frequency of the thin shell gravastar does not coincide with that of the 
black hole even when the radius of the thin shell gravastar is nearly its Schwarzschild 
radius. This is because different boundary conditions at the nearly gravitational radius 
are required for the quasinormal modes of the gravastar and black hole, (ii) The thin 
shell gravastar may be distinguishable from the black hole by using the quasinormal 
frequencies, (iii) The particular models of the thin shell gravastars are unstable against 
nonradial perturbations for the case of less compact models. 

The other of such potential data obtaining from the gravitational wave observatories 
is the values of the lowest few multipole moments of the black hole-like object (its 
gravitational mass M, its spin angular momentum J, its mass quadrupole moment 
Q) [9]. This data will be obtained through the so-called extreme mass ratio inspiral 
(EMRI) associated with the supermassive black hole candidate in the center of a galaxy. 
Only when the central black hole-like object is non-spherically symmetric, one may 
obtain non-zero values of J and Q. The quadrupole moment Q includes information 
about the matter properties and is determined hj Q = /M for the genuine black hole. 
Thus, one may, in principle, distinguish the genuine black hole from the alternative 
models to the back hole by using the values of M, J, and Q. 

So far, studies on the gravastar have mainly done for the spherically symmetric 
nonrotating case. As for the rotating gravastar, the first order rotational effect in 
the spin angular momentum J, which appears as the frame dragging effect, has been 
considered to investigate the so-called ergoregion instability of the gravastar mm- 
Since most black hole candidates are expected to rotate rapidly because of conservation 
of angular momentum and/or matter accretion, it is necessary to extend the study of 
rotational effects on the gravastar structure to the higher order case. In this study, 
we will therefore study rotating thin shell gravastars. As a first step, we focus on the 
case of slow rotation. We treat the rotational effects on the stellar structure as small 
perturbations about nonrotating spherically symmetric thin shell gravastar by using a 
method employed for investigation of the slowly rotating regular black hole m, which 
is based on the formalism to treat the slowly rotating star and the matching 

of the distinct spacetimes dSlES]. The rotational effects are taken into account up to 
the second order in the spin angular moment J. Including the second order rotational 
effects, we may estimate the value of the mass quadrupole moment Q of the gravastar. 
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which is a quite important quantity to distinguish the black hole from the gravastar 
as mentioned before. For the nonrotating unperturbed solution, we use the thin shell 
gravastar model considered by Visser and Wiltshire [3]. 

The plan of this paper is the following: In section 2 we show the basic setup for 
matching of the two distinct spacetimes. The master equations for slowly rotating 
thin shell gravastars are given in section 3. In section 4, we show numerical results of 
the rotational corrections to the structure of the thin shell gravastar. For the present 
gravastar models, the limit of the slowly rotating thin shell in vacuum asymptotically flat 
spacetime and the properties of the ergoreion are discussed in section 5. The conclusions 
are given in section 6. In this paper, we use the units of G = c = 1, where G and c are 
the gravitation constant and the speed of light, respectively. 

2. Formulation: Matching of the two distinct spacetimes 

2.1. Spacetimes glued with a thin shell 

To obtain rotating models of the thin shell gravastar, in this study, we match two 
distinct spacetimes appropriate for them with a timelike thin shell. Following Visser 
and Wiltshire [3], in the zero-rotation limit, the interior and exterior spacetimes of the 
thin shell are assumed to be the de Sitter and Schwarzschild spacetimes, respectively. 
We may obtain solutions of the nonrotating thin shell gravastar if the two spacetimes 
are properly matched at the coordinate radius r = R, where 2M < R < L with 2M 
and L being the Schwarzschild and de Sitter horizon radii, respectively. Then, the thin 
shell gravastars, in general, have a massive spherical thin shell aX r = R. 

In this study, we shall restrict ourselves to the case of slow rotation, i.e., e = 
<C 1, where 12 is the angular velocity of the thin shell measured by an observer at 
infinity and is the Kepler frequency of the nonrotating thin shell gravastar. Following 
a standard technique developed for the analysis of the slowly rotating star naimiii], 
we treat rotational effects on the stellar structure as stationary and axisymmetric 
perturbations around a zero-rotation solution. The rotational effects are taken into 
account up to an accuracy of O(e^) and we further assume that rotating thin shell 
gravastars have equatorial symmetry and time-azimuth reflection symmetry. Thus, the 
line elements for the interior and exterior spacetimes of the thin shell may be written 
in the common form, given by 

ds^ = dag dx°‘dx^ 

= -/(r)(l + 2Air, e))de + ' (l + dr 

f{r) V rf{r) ) 

R r^(l -I- 2e^k{r, 9)) [dO"^ + sin^ 9 {dcj) — euj{r)dt}'^] 0{e^), (1) 

where g^g stands for the component of the metric tensor. Here and henceforth, the greek 
letter subscripts and superscripts denote indices running from 0 to 3 for four-dimensional 




5 


Slowly rotating thin shell gravastars 


tensor quantities. Outside and inside the thin shell, respectively, we assume that 
(x+)/^ = (t+,r+,0+,0+), 

/(r) = /+(r+) = 1 - 

y^-t- 

(h(r, 6*), /c(r, 9),m{r, 9),u){r)) = 

(h+(r+, 6>+), A;+(r+, 6^+), m+(r+, 6^+), a;+(r+)), (2) 

and 

= (f,r",0“,0“), 

(r~)^ 

/W = /-(r-) = l-^, 

(h(r, 6^), /c(r, 9),m{r, 9),u){r)) = 

{h~{r~, 9~), k~{r~, 9~),m~{r~, 9~),uj~{r~)). (3) 

Here and henceforth, the quantities with the superscripts (■*■) and (“) mean the 
quantities dehned outside and inside the thin shell, respectively. For the sake of 
simplicity, however, the superscripts (■*■) and (“) are sometimes omitted when it does 
not cause any confusion. We further assume that 


h{r,9) = hoir) + h2{r)P2{cos9), 
m{r,9) = mo{r) + m2{r)P2{cos9), 

k{r,9) = k2{r)P2{cos9). (4) 


For the exterior spacetime of the thin shell, the functions cu’*', h'^, k~^, and m~^ are 
the solutions regular at inhnity obtained by solving the vacuum Einstein equations, see, 
e.g., mm, given by 


u 


+ 


2J 


m+ = 
h+ = 

ht = 


k+ 


6M - 




6M 


+ 




r — 2M r^(r —2M)’ 


r 


Z! 

„4 


Mr3 
B 


+ ^ + bq 


2M 




Ql 


r 

M 



m 2 = {r 


2M) 






(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 

( 10 ) 


where 6M means change in the total gravitational mass due to the rotation, J the total 
angular momentum, Q™ the associated Legendre function of the second kind, and B 
an integral constant. The integral constant B is related to the quadrupole moment of 
the gravastar Q through the relation Q = J^/M + 8BM^/5 (see, e.g., [I2])- Note that 
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the case of i? = 0 corresponds to the slowly rotating Kerr solution. 
Legendre functions and are, respectively, defined by 


2 / ^ ^(5 


Q 2 

Q^{x 




X) = 


+ 


3{x‘^ 


x^ 


1) , X + 1 
— log 


= 


x^ 


3x^ 3x , 

+ — log 

— 1 2 X 


X — 1 ’ 

X + 1 


The associated 


( 11 ) 

( 12 ) 


For the interior spacetime of the thin shell, the functions uj~, k~, and m~ are 
the solutions regular on the symmetry axis obtained by solving the Einstein equation 
with a positive cosmological constant, see, e.g., m, given by 


u — Cl, 

tuq = 0 , 

ho = C 2 , 

_ C 3 /-3L2 + 5r2 

2 “ 8r2 V L2/-(r) 

Cs f3L^ + 4r^ 
^ ~ 8r^L V L 
= -rf~{r)h2, 


3Lf~ (r) Arctanh(r / L)\ 
r ) ’ 

3(L^ + r^)Arctanh(r/L) \ 
r / ’ 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


where Ci, C 2 , and C 3 are integral constants. 


2.2. First junction conditions 


In order to match the two distinct spacetimes given in the previous subsection properly, 
the so-called Israel’s junction conditions [T51 [TB] have to be fulhlled. Israel’s junction 
conditions are composed of two parts. The first conditions are related to the induced 
metrics on the thin shell, h^ab, and the second conditions to the extrinsic curvatures on 
the thin shell, K^ab- Here, let us consider the first conditions, which require that h'^ab 
is continuous across the thin shell, i.e.. 


[[hab]] = 0 , 


(19) 


hab = Qage'je^ , e" 


dx°‘ 
dy°- ’ 


( 20 ) 


where e“ means the tangent vectors to the thin shell, S, and y°' denotes the intrinsic 
coordinate of the thin shell, whose equation is given in the form x“ = x“(?/“). Here 
and henceforth, the double square brackets [[g]] imply the jump in the quantity q on 
the thin shell, i.e., [[g]] = gl^ — gl^, and the roman letter subscripts and superscripts 
denote indices for the three-dimensional tensor quantities defined on the thin shell. In 
this study, the intrinsic coordinate of the shell is written by y°‘ = (T, 0, $) and the 
equations of the thin shell in the interior and exterior spacetimes of the thin shell are 
assumed to be 


(X^Y = (A^T, R + £^*(0), 0 + £^(/*)®(e), 4) + 0(Y), 


( 21 ) 














Slowly rotating thin shell gravastars 


7 


where R = is the coordinate radius of the thin shell in the zero-rotation limit. In 
this study, we may further assume that = 1 , A~ = constant, and (A)® = 0 by using 
the freedom in the coordinate choice (see m)- The tangent vectors to S may be then 
given by 

(e^)± = (7l±,0,0,0)+0(e3), 

(e^)^ = (0, 1 + e^{l^)%, O) + O(e^), 

(e^)± = (0,0, 0,1)+ 0(63). (22) 

To achieve separation of variables, we assume that 

e( 0 ) = eo + e 2 P 2 (cos 0 ), (23) 


where .^o and ^2 are constants. 

The e°- and e-c 
respectively, lead to 


The e°- and e-order first junction conditions, hni[[hTT]] = 0 and hrnelg))^^^]] = 0, 




riR) 


UJ- = Cl = 


2J 


f-{Ry " R^A- 

The conditions [[hee]] = 0 and [[h$$]] = 0 lead to 


iKoii - m-)% = 0. iKoii=0. 


^ + k2{R) 


= 0 . 


(24) 

(25) 

and the regularity 


We therefore have (/ )® = 0 because of the first two equations in 

on the symmetry axis. From the e^-order conditions of lim -5^ [[h^r]] = 0, we obtain 

£—>■0 2 


[[ho(i?)]] + 

[[h2{R)]] + 




+ 


Mto 


L^f-{R) R^f+{R) 

R^2 ^ 


= 0 , 

= 0 . 


(26) 


L^f-{R) R^f+{R) 

The conditions [[h-Te]] = 0 and [[h©$]] = 0 are automatically satisfied because they are 
smaller than e^-order quantities. Nonzero components of the induced metric may be, in 
terms of quantities associated with the exterior spacetime of the thin shell, given by 

hrr = - /+ - ^in^ 0(a;+)2 + (In + 2h+ 

+ {(In f~^y^2 + } P 2 ] + 0{ey , (27) 

hee = R^ + e‘^2R [^ 0 + + {Rkt + ^t)P2] + 0(e"), (28) 

= sin^ 0 / 1 ©© , (29) 

hT<s> = = —eR^oj^ sin^ 0 + 0(6^), (30) 


where q' means dq/dr\r=R. 

In summary, the first junction conditions give us the six constraint conditions, the 
six conditions given in equations fl2411 - fl26l) . for the eleven constants J, (5M, B, Ci, C 2 , 
C 3 , A~, and Thus, the first junction condition, ([19]), can be fully satisfied as 
shown later. 
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2.3. Second junction conditions 


While the first junction conditions, (|T^ . ensure that the two distinct spacetimes are 
joined smoothly, the second junction conditions give us information about the matter 
distribution on the matching surface S in terms of the jump in the extrinsic curvature 
across S (see, e.g., mm)- According to Israel’s junction conditions, the stress-energy 
tensor of the thin shell Sab is given by 


S.I, = ^ (l[-fCJ| - h.^[[K ]\), 




ol 0 


and K = habK 


ab 


(31) 

(32) 


where is the unit normal vector to S, for which n,Q,e" = 0 and = 1, and the 

semicolon (;) denotes the covariant derivative associated with gag. Nonzero components 
of are given by 

1 


Ur = 


V7 


^ 2 ^0 -F 7712P 2 (cos 6^) ^^^ 4 ^ 


rV7 


6 


He = -e^^dgP2{cose) 0{ep. 
V / 


(33) 

(34) 


Nonzero components of the extrinsic curvature up to an accuracy of 0(e) are then given 
by 

r 


Kj. = 


Plq — — 


2^7 
$ /7 


+ 0 ( 6 ^), 

+ 0(e^). 


R 


$ _ /(2a; -h Ru') - Ruf . , 3 ^ 


A/ = e 


2i?x/7 


-A + 0 (e^), 


R^u' 




K = 


2A^ 

^f + Rf 


O(e^). 


(35) 

(36) 

(37) 

(38) 

(39) 


2A/7 

Explicit forms of the e^-order coefficients of the extrinsic curvature component are 
summarized in Appendix A. 

The matter three-velocity tangent to S, and the total energy density measured 
by an observer with n“, a, are, respectively, defined by the timelike eigenvector and the 
corresponding eigenvalue of S^, given by 

u°‘Ua = —1, (40) 

where is related to the matter four-velocity associated with the thin shell through 
n" = The stress tensor of the thin shell, -jab, measured by an observer with is 

defined by 

'lab = S'^'^qacqbd, ( 41 ) 
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where qab is the projection tensor orthogonal to the dehned by qab = hab + UaUb- For 
a perfect fluid, ■jab is proportional to qab, lab = PQab, and p is interpreted as the isotropic 
pressure of the fluid. 


3. Master equations for slowly rotating thin shell gravastar 


3.1. Non-rotating thin shell gravastars 


In the limit of e ^ 0, fHOl) may be solved if we assume that = 




5^ + 0(e) with 


(5^ being the Kronecker delta. Thus, we obtain the energy density and pressure in the 
zero-rotation limit as follows: 

vr-vr 


an = lim a = 


Pq = limp = 
£->■0 


47ri? 



(42) 

(43) 


Note that the relation jab = PQab + 0{e‘^) may be trivially confirmed in the present 
situation. In order to specify an equilibrium solution of the nonrotating thin shell 
gravastar, following [3], we give values of two dimesionless parameters R/M and M/L 
under the condition of 2M < R < L. Then, we may obtain nondimensional fundamental 
physical quantities for the nonrotating thin shell gravastar, A~, Ana^R^/M, and Po/<^ 0 ) 
by using the first equation in (IM)) . (I4^ and (l43)l . 


3.2. e-order rotational effects: frame dragging 


In this study, we focus on the case of the uniformly rotating thin shell around the 
symmetry axis. Thus, is given by 

u^d^ p. 

= — = VL = constant, = 0 , (44) 

dT 

where hi is the angular velocity of the thin shell and = 0(e). As a result, the 0 
component of (140|) is automatically satished since Af = O(e^) = S'®, and the T and $ 
components of (HHl) become 

+ Sin = -a + O(e^), (45) 

and 

^| + ^|f] = -aof^ + 0(e3), (46) 

respectively. Equation fl46D leads to the relation between the angular velocity of the 
thin shell hi and the total angular moment J, given by 


= eflfc 


+ CTO 

J /r+27F 
T" TFfl - (fl - 3M)yF ’ 


( 47 ) 
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or 


J 


.ZL \ U If • 

\/F + 2x/F 


(48) 


Here, as mentioned before, is the Kepler frequency of the nonrotating thin shell 
gravastar, dehned by y^M/R^. Equation fl45|) determines the total energy density up 
to an accuracy of O(e^). As seen from the above, we conhrm that given in (14^ is a 
timelike eigenvector of within an accuracy up to O(e^). The angular velocity of the 
frame dragging in the interior spacetime of the thin shell u~ is, in terms of J, given by 
the second equation in fl241) . 


3.3. -order rotational effects: deformation of the matter distribution due to rotation 


As shown in the previous subsection, the three-velocity of the thin shell is given 
by fl44D . Then, the corresponding eigenvalue a is given by fl45|) . a is the total energy 
density and may be written in the form, given by 

a = ao €^5a + 0 {e^) 

= (To -h e^((5cro -h (5(72^2) + 0{e ^), (49) 

where the explicit forms of 6 ao and 6 a 2 are given in Appendix B. While H is related to 
J through (|47|) . as mentioned before, u'^ is determined by the normalized condition of 

habU’^u'^ = -1 . (50) 


Thus, we have 


T 

u = 


u = 




+ e^ 


'-(/+)'^+(0) + R^iyik - 0J+) sin^ 0 h+(i?, 0)' 


2/7+' 






0 {e^) , 
H + 0(e=^) 


(51) 

(52) 


By using these expressions for we may obtain nonzero components of the stress 
tensor of the thin shell, 7 ^, given by 

7t =Pogr + C»(e^)) 7$ =Pog$C>(e'), 1% = Poq% + 0{e^), (53) 

r 2(flfci?3 - 2 J) / JL2 - J(3M - 2R) + MR^Vtk 


7 =Po- 


87ri?2 \ 3i^3/^ 




+ 




- )' - 


(3M2 - 3Mi? + R^)if + R{R - M)m^ 


e^P2 {2{JL^ - R^nk){R^nk-2J^ 


SnR^ {' 


3Lm^^f^ 


+ R^[[^/fiK + K)]] 


(3i?2 - 2L2)^2" + - ‘2R^)mf ^ 2(2J - R^nk){J{3M - 2R) + MR^Vtk) 


L^y/T 

(3M^ + 3MR - 2R^)^+ + R{R - M)mt 


3R^^/J+ 


+ 0{e% 


(54) 
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7 


e^sin^^ / 2(2 J - 

167r/?2 I R^f+ 
J(3M - 2i?) + 

WF y 


/ JL‘^-R^Vtk 

V iF/r~ 



0 {e% 


(55) 


where 7 =^ = ( 7 I ± 7 |)/ 2 . 

In this study, we focus on the thin shell composed of a perfect fluid. For a perfect 
fluid, as mentioned before, the stress tensor satishes the relation 7^6 = p qah- In equations 
fl5^ - (l5^ . we see that the relation 7^6 = pQab is satished within an accuracy up to O(e^) 
except for 70 and 7 I,. If the thin shell matter is described as a perfect fluid, 7 I, and 
70 are written as 

$ $ 

= (po + e^<5p)(l+ «'*’«$) +O(e^) (56) 

= Po + e^5p + PqU^u^ + 0{e'^), 


and 


1 % ^pq% ^Po + e‘Sp + 0(e'‘). 


(57) 


Thus, the condition for the thin shell to be composed of a perfect fluid is given by 


$ 

7$ 


e $ 

7© = PoM Uq>. 


This condition may be rewritten in the form, given by 


[[F. 


/4J2 2J{M + R)nk 


+ MR^Q^ 


3//T [r3 r 

V 3 / 


(58) 


(59) 


The pressure of the thin shell p may be given by 

p = 7 + - 1 pan*Ui 
= Po + e^3p 

= Po + €^(5po + ips-Pe), 


(60) 

(61) 

(62) 


where the explicit expressions for (5po and 5 p2 are summarized in Appendix B. 

Since the perfect-fluid thin shell is considered, the total energy and pressure 
perturbations have to be related by the relationship, 

6 p = ^ Sa , (63) 


dp 

where — is a function of the thermodynamic quantities. This implies that 
da 


X - X 
opo — —— oao : 

da 
dp . 


dp2 = ^da2. 

da 


(64) 


Once the coefficient is known, we have two more equations for determining the eight 
da 

constants 5M, B, C 2 , O 3 , and associated with e^-order rotational effects. 
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The second equation in flMD . dSU]) . the second equation in (12^ . and the third 
equation in fl25|) are composed of four coupled linear algebraic equations among the four 
constants B, C 3 , and associated with e^-order quadrupole perturbations. Thus, we 
may obtain solutions of the quadrupole perturbations by solving these coupled linear 
algebraic equations (see Appendix C). 

To determine the four constants 5M, and C 2 associated with the spherically 
symmetric perturbation, an additional condition is required because we have so far 
obtained the three algebraic equations; the hrst equation in (IMl) . the hrst equation in 
(|26|) . and the second equation in ([25]). As a fourth condition, in this study, we impose 
that 6ao = 0. By adding this condition, we may uniquely determine the four constants 
SM, and C 2 . 

From the hrst law of thermodynamics, changes in the particle number density of 
the thin shell 6n is given by 

6 n = — — —6 a, (65) 

(^o+Po 

where uq is the number density in the zero-rotation limit. Changes in the total particle 
number of the thin shell 6 N is given by 

6N = no [ ( - + 6u^ + v ?"^ \/—h dQd^, (66) 

J + Po \/—h ) 

where h means the determinant of hah- This leads to 

N R 3R^f+ ’ ^ ^ 

where N denotes the total particle number of the nonrotating thin shell gravastar. Note 
that 6 N is given only by the spherically symmetric perturbation quantities, because the 
terms related to the quadrupole perturbations vanish after the angular integration. 

4. Numerical results 

4.I. Thin shell gravastar in the zero-rotation limit 

Mazur and Mottola have proposed the original gravastar model, in which a hnite- 
thickness shell described by the stiff equation of state, p = a, is assumed to be around 
the radius at which the horizon would have formed in the standard classical general 
relativity [2] . A simplihed version of the gravastar model has been considered by Visser 
and Wiltshire [3]. This simplihed model is called the thin shell gravastar because the 
hnite-thickness shell of the original model is replaced by the single thin shell. They 
have made detailed investigations of the properties of the thin shell gravastar and shown 
that there are gravastar models stable against spherically symmetric perturbations. If 
spherically symmetric nonrotating models are unstable, to study rotational effects on 
such models will lose its physical signihcance. In this study, therefore, we only consider 
the stable thin shell gravastar models as solutions in the zero-rotation limit. The stable 
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Figure 1. The energy density of the thin shell ctq in the zero-rotation limit as a 
function of R/M. 



Figure 2. The same as figure [T] but for the ratio of the pressure of the shell to the 
energy density of the thin shell po/ao- 
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Figure 3. The pressure of the thin shell po in the zero rotation limit as a function of 
the energy density of the thin shell CTo- 
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Figure 4. The adiabatic index T for the thin shell matter as a function of the energy 
density of the shell cto. 


nonrotating thin shell gravastar models that we employ in this study correspond to the 
case of 1/ = 0 in [3] . 

As mentioned before, we assume that 2M < R < L to avoid appearing solutions 
with horizons. Following Visser and Wiltshire [3], we require that the thin shell matter 
satisfies the weak and dominant energy conditions. The weak energy condition uo > 0 
and 02] lead to the constraint equation, given by 


M - V2M2 J 


( 68 ) 


Thus, we have 


R 

2 < T7 < 

M “ 



(69) 


Although it is a complicated numerical task to translate the dominant energy condition 
for the thin shell matter, (Tq > |po|) fo fhe condition given explicitly in terms of R, M, 
and L, a simultaneous analysis of the dominant energy condition and the radial stability 
condition can derive the condition, 

— < 0.024304, (70) 

for the thin shell to satisfy the both conditions (see, for detailed derivation, 0). In this 
study, therefore, we consider solutions of the thin shell gravastar in the zero-rotation 
limit characterized by fl6^ and flTD]) . 

To investigate rotational effects on the structure of the thin shell gravastar, changing 
values of i?/M, we obtain the five sequences of equilibrium solutions of the spherical 
gravastar characterized by the fixed values of M/L, which are chosen as the following: 
M/L = l/\/^, 1/5, 1/6, l/\/40, and 1/7. For showing numerical results of the 
gravastar, in this paper, we will use the units of L = 1. In figures [D and E] the total 
energy density Uo and the ratio of the pressure to the total energy density Po/aQ of 
the thin shell for the equilibrium sequences of the gravastar are given as functions of 
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R/M, respectively. In these figures, each line corresponds to an equilibrium sequence 
characterized by the fixed values oi M/L whose values are given in these figures. In 
these figures, we see that the energy density of the shell do is a decreasing function of 
R/M. As can be seen from figure [21 the endpoints of each sequence considered in this 
study correspond to the solutions whose thin shells are described by the stiff equation 
of state, po = <7o- For a sequence of equilibrium solutions characterized by a fixed value 
of M/L, therefore, there are two solutions with the stiff thin shell. Although in all the 
solutions shown in figures [H and [21 the thin shells satisfy the weak and dominant energy 
conditions and are stable against radial perturbations, i.e., they are equally plausible, 
thin shell gravastar solutions having smaller radius would be more relevant because they 
should be black hole-like objects. Note that most unperturbed solutions we consider 
have their radius of R/M < 3. Such solutions might be nonlinearly unstable as pointed 
out in [m HE] . 

Considering the equilibrium sequences characterized by fixed values of M/L, we 
may obtain equations of state for the thin shell matter, po = Lo(c’‘o), depending on 

values of M/L. These equations of state are shown in figure [3l By assuming that these 

dp 

equations of state specify values of — appearing in f[63D . we may reasonably determine 

dcT 

the equation of state for the perturbation. We therefore assume that 
da 


duo 

OR 


(71) 


2(a-o -hpo) OR’ 

The adiabatic index T for the equations of state given in figure [3] is shown in figure [H 
Here, the adiabatic index T is defined by 

( 72 ) 

Po da 

Since we obtain the nonrotating solutions of the thin shell gravastar given in figures 
[U and [21 we may evaluate rotational effects on the structure of the spherical gravastars 
as small perturbations around them. In this paper, we shall basically show numerical 
results of the perturbation quantities for e = 1 or hi = 


4 -2. e-order rotational effects on the thin shell gravastar 

The e-order rotational effects appear as the frame dragging effects in the exterior 
and interior spacetimes of the thin shell as discussed before. An observer at infinity 
can measure the frame dragging effects in the exterior spacetime as the total angular 
momentum of the gravastar J. In figure [5l we show the dimensionless total angular 
momentum of the thin shell gravastar J/M"^ along the equilibrium sequences as a 
function of R/M. In this figure, we observe that the dimensionless angular momentum 
for the M/L constant equilibrium sequences is basically an increasing function of R/M, 
however, it becomes a decreasing one for some equilibrium sequences when the gravastar 

is very compact (when R/M < 2.25). In figure [6l values of —-—- for the equilibrium 

uj~^{R) 

sequences given in figures [T] and [21 are given as functions of R/M. 
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Figure 5. The dimensionless angular momentum of the gravastar, as a function 

of R/M. 



R/M 


Figure 6. The same as figure [5] but for 

u}+[R) 


4-3. -order rotational effects on the thin shell gravastar 

First, we show physical quantities associated with quadrupole perturbations. Of 
these quantities, the quadruple moment is of the largest signihcance, because it 
represents information about the matter distribution that does not appear for spherically 
symmetric solutions and can be measured by analyzing motions of a test particle 
around it. As argued before, the slowly rotating Kerr black hole is characterized by 
the particular relation of the quadrupole moment QM/J"^ = 1. As a result, we may 
conclude that slowly rotating compact objects are not black holes if their quadrupole 
moment does not obey the relation QM/J‘^ = 1. 

In hgure [71 we show the inverse of the dimensionless quadrupole moment 
{QM /and the dimensionless quadrupole moment QM /for the hve equilibrium 
sequences characterized by a hxed value of M (see section 4.1). In the left panel of this 
hgure, we observe that values of {QM /vanish at a single point along the three 
equilibrium sequences with M = 1/6, l/\/40, and l/\/7. This means that because 















Slowly rotating thin shell gravastars 


17 



Figure 7. Left: the inverse of the dimensionless quadrupole moment, {QM/J^)~^, 
as a function of R/M. Right: the dimensionless quadrupole moment, QM/J"^, as a 
function of R/M for the range of 2.1 < R/M <2.8. 



Figure 8. The square of the ellipticity of the thin shell, e^/e^, as a function of R/M. 



Figure 9. The same as figure [8] but for the normalized quadrupole component of the 
radial displacement of the thin shell, / R. 
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R/M 

Figure 10. The same as figure [8] but for the normalized quadrupole component of the 
radial displacement of the thin shell, 


perturbation quantities diverge for these unperturbed solutions, the nonrotating thin 
shell gravastar solutions with {QM/J^)~^ = 0 cannot rotate with the perturbative 
approximation under the assumptions introduced in this study, i.e., assumptions of the 
perfect-fluid thin shell, fl58l) . and the equation of state for the thin shell, flTTl) . The 
same singular behavior occurs in all the other quantities associated with quadrupole 
perturbations for the unperturbed solutions with {QM/.P)~^ = 0. Existence of 
these solutions with {QM/.P)~^ = 0 necessarily implies that the present perturbative 
approach also becomes irrelevant around such solutions. As seen from the left panel 
of figure [71 however, nonrotating solutions with sufficiently lower and higher values of 
R/M than those of the ones that cannot rotate with the perturbation approximation do 
not show any difficulty for investigation of rotational effects on their structure with the 
present treatment. Since we are interested in black hole like objects, hereafter, we will 
focus on the cases of the nonrotating models with R/M < 2.8, for which no singular 
behavior appears in the perturbation solutions. 

In the right panel of hgure [71 we show the dimensionless quadrupole momeut 
QM/,P for the hve equilibrium sequences characterized by a hxed value of M and 
R/M < 2.8. From the right panel of hgure [71 we can see that some solutions have 
the same dimensionless quadrupole moment as that of slowly rotating Kerr black holes, 
i.e., QM/J'^ = 1. For these solutions, we cannot distinguish slowly rotating gravastars 
from slowly rotating black holes by using the quadrupole moment. Another remarkable 
feature of the quadrupole moment of the slowly rotating thin shell gravastar is that for 
sufficiently compact models, values of QM/J'^ are less than unity, which are usually 
attributed to the matter distribution deformed prolately. For standard self-gravitating 
rotating objects, their matter distribution is oblately deformed due to the centrifugal 
force and their values of QM /are larger than unity. In order to see degree of 
the deformation of the thin shell associated with its shape, in hgure [H we show the 
normalized square of the ellipticity of the thin shell e^/e^ for the hve equilibrium 
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R/M 


Figure 11. The same as figure [8] but for the normalized quadrupole component of the 
energy density perturbation of the thin shell, 802 !<Jq. 


sequences characterized by a fixed value of M, where e^/e^ is defined by 




kt[K) + ^ 


(73) 


Note that we have 


HR) +1 


0. With this definition, positive (negative) values 


of e^ mean that slowly rotating thin shells are of oblate (prolate) shape. From figurelHl we 
see that for the case of R/M < 2.8, the solutions belonging to the equilibrium sequences 
with M = 1/\/M and 1/5 always have positive values of the square of the ellipticity of 
the shell, while for the solutions belonging to the other equilibrium sequences, values of 
e^ are basically negative but becomes positive at lower and higher values of R/M. Thus, 
some slowly rotating thin shell gravastars have an oblate shape thin shell but their values 
of QM/.P are less than unity. In figures [9] and [101 we show the dimensionless radial 
displacements of the spherical shell and ^2 /R, respectively. For the nonrotating 

gravastar models with 2.1 < R/M < 2.8, we can see that values of ^2 /^ ^ 2 /^ 

are always positive. The energy density and pressure perturbations rescaled by their 
unperturbed quantities, ^ 02 / 0 ^ and hp 2 /po, are given as functions of R/M in figures fTTl 
andHH respectively. We observe that values of 5 a 2 /o'Q and 8 ^ 2 /Vo are always negative 
in the range of 2.1 < R/M < 2.8. Values of Cs are shown in figure [l3] along the 
five sequences of the nonrotating thin shell gravastars. From this figure, it is found 
that Cs < 30 for all the nonrotating solutions for 2.1 < R/M < 2.8. Since all the 
perturbation quantities given before are of the order of 0(1), we can conclude that our 
perturbation approach is suitable for the unperturbed solutions with 2.1 < R/M < 2.8. 

Next, we show the physical quantities associated with spherically symmetric 
perturbations. Since Sag = 0 is assumed (see section 3.3) and values of dp/da are 
finite as shown in figure [3l we have 5pQ = 0. In figures [H] and [HI respectively, we 
show relative changes in the total gravitational mass, 5M/M, and in the total particle 
number, 8N/N, of the thin shell gravastar due to rotation for the five equilibrium 
sequences characterized by a fixed value of M. In these figures, it is found that their 
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Figure 12. The same as figure[8]but for the normalized quadrupole component of the 
pressure perturbation of the thin shell, ^vilvvs- 



Figure 13. The same as figure|S]but for C 2 - 


behaviors are qualitatively similar to each other, that their values are always positive, 
and that there are the maximum values of 5M/M and 5N/N for the four equilibrium 
sequences with M = 1/7, l/\/40, 1/6, and 1/5. The dimensionless radial displacement 
^Q/R = ^^/Ris shown in hgure [T6l From this hgure, we see that ^o/R's as functions of 
R/M show similar properties to those of 5M/M and 6 N/N. Values of C 2 as functions 
oi R/M along the five equilibrium sequences are shown in figure fTTl In this figure, we 
see that values of C 2 are negative and of the order of 0(1). 

5. Discussions 

5.1. Slowly rotating thin shells with isotropic pressure: limit of L ^ 00 

In section 3, we give master equations for obtaining slowly rotating thin shell gravastars. 
This formalism can also be applicable in the case oi L ^ 00 . In this limit, the de Sitter 
core vanishes and the place occupied by a de Sitter core becomes vacuum. In other 
words, only a slowly rotating thin shell exists in a vacuum and the thin shell matter is 
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Figure 14. The same as figure |8] but for the relative change of the gravitational mass 
6 M/M. 



Figure 15. The same as figure [5] but for the relative change of the total particle 
number of the thin shell, SN/N. 



Figure 16. The same as figure [8] but for the normalized spherically symmetric 
component of the radial displacement of the thin shell ^o/R- 
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Figure 17. The same as figure |8] but for C 2 - 



Figure 18. Energy density ctq and pressure po of the thin shell in vacuum 
asymptotically flat spacetime in the zero rotation limit as functions of R/M. 


the only source of the gravity. This is another interesting and fundamental situation 
apart from the gravastar. Therefore, similar but not the same situations have been so 
far considered by several authors, e.g., liniisDiEiiEg. 

The quadrupole metric perturbations for the interior spacetime of the thin shell, 
given in fllbp and fll7p . vanish in the limit of L —)■ oo. However, this is not appropriate 
behavior for the metric perturbations inside the slowly rotating thin shell, because in 
general, the quadrupole metric perturbations need not vanish there. In order to obtain 
physically suitable solutions in this limit, we will use another integral constant Cs, 
dehned by Cs = 5L^ C 3 . In the limit of L —)■ 00 , we then have appropriate solutions, 
given by The integral constant C 3 is determined with the same 

procedure described in the previous sections. For all the other perturbation quantities, 
we may take the L ^ 00 limit without any difficulty. 

Since the limit of L —)■ cxd is taken, it will be convenient to use the physical quantities 
nondimensionalized by a length scale of M. Henceforth, we therefore employ the units 
of M = 1 to show the numerical results for the slowly rotating thin shell. For the weak 
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Figure 19. The same as figure[T8]but for the dimensionless angular momentum J/M^ 
and u)~/uj^{R). 




Figure 20. Left: the dimensionless quadrupole moment of the thin shell, QM/ J^, as 
a function of R/M. Right: the quadrupole components of the perturbation quantities, 
^2 1^1 ^2 1^^ 5p2/po, ^cr 2 /cro, e^/e^, and C 3 as functions of R/M. 




Figure 21. Left: the spherically symmetric components of the perturbation quantities, 
^o/R, SM/M, and SN/N as functions of R/M. Right: C 2 as functions of R/M. 
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energy condition, (Tq > 0, to be satisfied, for the nonrotating thin shell, we require the 
inequality [68l In the limit of L —)■ cx), this condition becomes trivial and the thin shells 
with (To > 0 are obtained for R > 2M. In this paper, however, we focus on the cases 
of i? < lOM for the sake of conciseness. In figure [HI the energy density and pressure 
of the nonrotating thin shell are shown as functions of the coordinate radius of the thin 
shell R. These solutions are unperturbed ones for investigating rotational effects on the 
structure of the spherical thin shell with the perturbation approach. From fiT2|) and 
(jl3]), we may see that the dominant energy condition ao > |po| for the nonrotating thin 
shell is met if 25M/12 < i? or 2.08M < i? is satisfied, which is consistent with results 
given in figure UHl In figure UHl the angular momentum of the thin shell J and the 
normalized angular velocity of the frame dragging inside the thin shell u~/u~^{R) are 
given as functions of R/M. The basic properties of those quantities are similar to those 
of the gravastar, as can be confirmed in figure [5l The dimensionless angular momentum 
J/M^ is basically an increasing function of the normalized radius of the shell R/M but 
becomes a decreasing function for sufficiently compact thin shell models. 

The physical quantities associated with quadrupole perturbations are summarized 
in figure 12(11 The dimensionless quadrupole moment QMj and the other perturbation 
quantities are shown as functions of R in the left and right panels of figure [201 
respectively. In the left panel of this figure, we observe that the dimensionless quadrupole 
moment QM/ is always smaller than unity and becomes negative when the radius of 
the thin shell becomes sufficiently large, i.e., when R > 2.9M. As can be seen in the right 
panel of figure [20l the normalized square of the ellipticity e^/e^ is basically negative, 
which means the shape of the slowly rotating thin shell is prolate. The dimensionless 
quadrupole radial displacement of the thin shell /R is always positive, while the 
quadrupole energy density and pressure perturbations are always negative. 

In figure [211 we show physical quantities of spherically symmetric perturbations 
6 M/M, 6N/N, ^o/R, (in the left panel) and C 2 (in the right panel) as functions of 
R/M. In this figure, we observe that the relative changes in the total gravitational 
mass 5M/M and the total particle number 6 N/N show similar behavior each other 
and are always positive. The dimensionless radial displacement of the thin shell ^o/R 
takes positive values regardless of the nonrotating thin shell models. In the right panel 
of figure EH it is found that the integral constant C 2 , which determines magnitude of 
metric perturbations inside the thin shell, increases rapidly as the radius of the thin shell 
R/M increases. This means that with perturbation approaches, rotational effects on the 
structure of the thin shell cannot be investigated if the radius of the thin shell becomes 
sufficiently large because for such thin shells with sufficiently large radii, infinitesimally 
small angular velocity leads to a thin shell rotating with finite rotation velocity. 

For a slowly rotating thin shell with isotropic pressure, an interesting limiting case 
is the Newtonian limit, which is obtained in the limit of M/R —)■ 0. In the Newtonian 
limit, we have 


M 


( 74 ) 
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(76) 
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(77) 
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2R 


(78) 


These results in the Newtonian limit are consistent with numerical results given in 
figures fT8] - [20l 

With perturbation approaches, as mentioned before, self-gravitating slowly rotating 
thin shells were studied by several authors, e.g., [191 1201 1211 122] • De la Cruz and 


Israel [20] and Pfister and Braun [21] obtained slowly rotating thin shells very similar 


to those discussed in this subsection. To investigate what is a possible source of the 
Kerr spacetime, de la Cruz and Israel obtained slowly rotating thin shells. Their 
main assumptions are the following: (i) the exterior spacetime of the thin shell 
is the Kerr one, (ii) the radial displacement in the exterior spacetime is given by 
= kR{a/— r+/R){l — r_/i?) cos^0, where fc is a constant, R is the radius 
of the nonrotating thin shell, a is the Kerr parameter, r_|_ and r_ are, respectively, the 
outer- and inner-horizon radii of the Kerr solution, and the radial coordinate r is of 
the Boyer-Lindquist coordinates, (iii) there is no change in the total gravitational mass 
due to rotation {6M = 0). Note that all the assumptions of de la Cruz and Israel’s 
are different from ours. Their solutions are uniquely determined in terms of the total 
gravitational mass M, the Kerr parameter a, the radius of the thin shell R, and the 
constant k, which determines Thus, both the oblate and prolate shape thin shells are 
equally possible for their solutions. Making the assumptions: (i) uniformly rotating thin 
shells (n = constant), (ii) the asymptotically fiat spacetime, (iii) the interior spacetime 
of the thin shell is fiat (Cs = 0), (iv) there is no change in the total gravitational mass 
due to rotation {SM = 0), (v) the radial displacement is given by cx sin^ 0, Pfister 
and Braun investigated rotational effects on the interior fiat spacetime of the thin shell 
and on the thin shell structures. Note that Pfister and Braun’s assumptions (iii)-(v) are 
different from ours. All the slowly rotating thin shells obtained by Pfister and Braun are 
of prolate shape. They also showed that even in the Newtonian limit {M/R —)■ 0), the 
shape of their slowly rotating thin shells is prolate. The common features of the slowly 
rotating thin shells obtained by De la Cruz and Israel 121 and Pfister and Braun IZH 
are that the thin shell matter stress is described by anisotropic pressure. 

As for the slowly rotating thin shells obtained in the present study, as argued before, 
the shape of the thin shells is basically prolate except for the sufficiently compact thin 
shell cases. Like in the case of Pfister and Braun [21], the slowly rotating thin shells 
obtained in this study have the prolate shape even in the Newtonian limit [see figure 
[20l (177|) . and (178]) ]. As shown so far, the slowly rotating self-gravitating thin shells 
frequently have prolate shape even in the Newtonian limit. Thus, we may guess that this 
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Figure 22. gtt on the constant 0 surfaces as functions of r/M for the model with 
M = 1/-\/2T, R = 2.244M, and = 0.4. Results for 9 = 0, 7r/4, and 7r/2 are given in 
the left top, right, and left bottom panels, respectively. The solid curves correspond 
to the results given by the line element (HD and the dashed curves to the results given 
by the line element (IMl) . 


occurrence of the slowly rotating thin shells with the prolate shape cannot be attributed 
to general relativistic effects like in the case of the centrifugal force reversal near the 
Schwarzschild black hole event horizon, see, e.g., [23]. It is intuitively-plausible that in 
the Newton mechanics, self-gravitating objects become oblate as the rotation angular 
velocity increases if only the centrifugal force is additionally exerted on the matter. In 
the construction of the slowly rotating thin shell, however, we employ the matching 
of the two distinct spacetimes. In this means, the distribution of the inhnitesimally 
thin matter is suitably determined for the two spacetimes to be glued continuously at 
a boundary hypersurface assumed particularly. Thus, it is not clear what is actually 
causing the occurrence of the prolate shape slowly rotating self-gravitating thin shells. 

5.2. Ergoregion of the slowly rotating thin shell gravastar within an accuracy up to the 
order of 

In general relativity, sufficiently compact and rapidly rotating objects can in principle 
have the so-called ergoregions, in which no physical observer with r = constant, 
6 = constant, and (f = constant can exist. Friedman has shown that an object having 
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Figure 23. The same as figure [22] but for the model with M = 1/7, R = 2.11M, and 
= 0.28. 

an ergoregion but no horizon in stationary asymptotic flat spacetime is unstable or 
marginally unstable to massless scalar and electromagnetic perturbations [21]. The 
same instability, of course, occurs for gravitational perturbations (see, e.g., 0 ). Since 
it is secular and basically has very long growth time, the ergoregion instability can 
substantially operate only when it has sufficiently short growth time. Thus, calculations 
of the growth time for the ergoregion instability in the system considered are required to 
see whether the ergoregion instability is driven effectively. However, general relativistic 
perturbations of fully stationary axisymmetric objects have in general been hard to deal 
with. Thus, a standard tractable way to treat ergoregions and ergoregion instability is to 
rely on an approximate method as follows: (i) Rotational effects are taken into account 
with a slow rotation approximation, (ii) Results obtained for slowly rotating objects 
are then extrapolated to the case of rapid rotation, (iii) The extrapolation results are 
employed for investigation of properties of the ergoregions and ergoregion instability. 
Schutz and Comins employed this approximation method for the first time to study the 
ergoregion and ergoregion instability of rotating stars [2^ [26] . The same method was 
applied for studies on the ergoregion instability of rotating gravastars mm- 

The ergoregion of stationary spacetime is defined by the region satisfying the 
inequality g^ut^t'-' > 0 with t^ being the time Killing vector. In the slow rotation 
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approximation used in this study, this condition leads to 

= gtt = -f{r) (l + 2e^h) + sin^ 6* > 0 . (79) 

If equation [79] is satished for some values of e satisfying 0 < e < 1, following Schutz 
and Comins [251 [26], we approximately regard that there is the ergoregion in the 
spacetime considered. In their calculations, Schutz and Comins further neglect the 
term 2e^h in [79] from physical considerations, although it is not strictly a consistent 
approximation [25l [26] . In most studies on the ergoregion instability of rotating 
objects, this approximation method devised by Schutz and Comins, in which h = 0 
is assumed in ([79l) . have been widely used, see, e.g., [TUTOllTT]. For the linear oscillation 
analysis, Schutz and Comins assumed that linear waves propagate through the spacetime 
described by the full metric, given by 

ds^ = —f{r)dt^ + g{r)dr‘^ + r'^dd'^ + sin^ 6{d4> — e 0dtY , (80) 


where / and g are the metric functions in the non-rotating limit. 

In the spacetime described by the line element o. there is a useful particular 
frame that is called a locally non-rotating frame (LNRF) by Bardeen, Press, and 
Teukolsky [28]. Observers associated with the LNRF follow trajectories dehned by 
d<j) = eoodt, dr = 0, and dO = 0. Thus, along trajectories of the LNRF observers, the line 
element is reduced to ds^NUF = “/(^) {1 + 2e^(ho(r’) -|- h2{r)P2{cos9))} dt^. Just inside 
the black hole event horizon, (Jslnrf becomes positive. This means that no physical 
observer with dr = 0 and d6 = 0 is allowed inside the event horizon and that descent of 
physical observers along r necessarily occurs there. Since we are concerned with objects 
with no event horizon nor spacetime singularity, in this study, we will consider the case 
of dsLNRF < 0 only. 

In the approximate treatment by Schutz and Comins, we have dsl^^^p = —f{r)dt^. 
Therefore, (Jslnrf < 0 is automatically satished regardless of values of e for objects with 
/(r) > 0. When Schutz and Comins’s approximate treatment is not adapted, or when 
the line element ([T]) is assumed, ds^NUF becomes positive for larger values of e in some 
gravastar models considered in this study. We observe that for such gravastar models, 
the region of (Jslnrf > 0 appears on the symmetry axis as values of e increase. 
Thus, the region of (Jslnrf > 0 appears if ho + h 2 < 0, and the condition for occurrence 
of dsLNFF > 0 is given by 


e2> 


2(ho + h2) 


if Hq /12 < 0 . 


(81) 


From this inequality, we may obtain the condition for no appearance of the dsp^^pp > 0 
region, given by 


e2< 




ii ho + h2 <0 , 


(82) 


Max]— 2{ho -\- ^ 2 )] 

where Max[g(r)] means the maximum value of the function q{r). In this study, we 
require this condition for values of (Jslnrf tie negative everywhere. For the present 
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Figure 24. Inner and outer radii of the ergoregion on the equatorial plane for the 
models with M = 1/y/^ and R = 2.244M. The two intersection points between 
the solid or dashed curve and a eJ/M"^ constant line are the radii of the ergoregion 
boundary. The solid curve corresponds to the result given by the line element © and 
the dashed curve to the result given by the line element 



gravastar models, we observe that the function + h^l has the maximum value at 
r = R. 

To see properties of the ergoregion of the gravastar, the two models are examined: 
the models having {M,R) = (l/\/^, 2.244M) and {M,R) = (1/7,2.11M). The 
former is the most compact model among the most massive equilibrium sequence (the 
M = l/\/2l equlinrium sequence), and the latter is the most compact model among 
the gravastar models calculated in this study. In hgures [22] and [23l hxing a value of 9, 
we exhibit g^s as functions of r/M for the gravastars with (M, R) = (l/\/^, 2.244M) 
and {M,R) = (1/7, 2.11M), respectively. For the models given in figures W2\ and l23l 
the values of are assumed to be = 0.4 < el 0.449 and = 0.28 < el 0.289, 
respectively. In these figures, the solid and dashed curves correspond to the results 
given by the line elements ([I]) and (|80|) . respectively. The left top, right and left bottom 
panels in figures [22] and [23] show the results for 6^ = 0, 7r/4, and 7r/2, respectively. Note 
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that in these hgures, gus> for r < R are different coordinate components from those for 
r > R, because the coordinate functions inside the thin shell are defined independently 
from the ones outside the thin shell, i.e., gu = gt-t- for r < i? and gu = gt+t+ for r > R. 
Thus, jumps in g^s appearing at r = i? are of spurious structure and have no physical 
meaning. We observe in these figures that values of gu decrease everywhere when the 
2e‘^h term in ([1]) is neglected. This implies that for the gravastar considered in this study, 
use of the line element flHOj) results in a smaller ergoregion. From figures [22] and [23l it 
is found that there is no ergoregion inside the thin shell. We also see from these hgures 
that values of gu are negative along the symmetry axis (i.e., there is no ergoregion on 
the symmetry axis) and that the cross-section of the ergoregion perpendicular to the 
symmetry axis is the maximum on the equatorial plane like in standard rotating objects. 
In hgures [24| and [25l we show the boundaries of ergoregion on the equatorial plane for 
the models having (M, i?) = (l/\/^, 2.244M) and {M,R) = (1/7, 2.11M), respectively. 
In these hgures, again, the solid and dashed curves correspond to the results given by 
the line elements ([1]) and fl80l) . respectively. The two intersection points between the 
solid or dashed curve and a eJ/M^ constant line are the radii of the ergoregion boundary 
for the gravastar rotating at a given value of eJ/M^. The minimum value of eJ/M^, 
^minJ/M ‘^1 for each curve indicates the value of the dimensionless angular momentum 
at which the ergoregion hrst occurs. The maximum values of eJ/M^ in hgures |24| and 
m are given by values of ecJ/M"^ determined by fl82l) . Note that gt-t-^s are always 
negative for the models considered in this study, so that the inner boundary of the 
ergoregion coincides with a radius of the thin shell. From hgures [2l| and [25l we see that 
the ergoregion forms at a smaller value of eJ/M^ when the line element ([T]) is assumed. 
The parcentage diherences between emin<7/Tf^’s given by the line elements ([I]) and fISOll 
are approximately 24.4% and 12.8% for the models with {M,R) = (l/\/2l, 2.244M) 
and {M,R) = (1/7, 2.11M), respectively. These relative diherences become smaller as 
the gravastar model becomes more compact as pointed out by Schutz and Comins |2S]- 

6. Conclusions 

We have constructed the solutions of slowly rotating gravastars with a thin shell. In 
the zero-rotation limit, we consider the gravastar composed of a de Sitter core, a 
thin shell, and Schwarzschild exterior spacetime. The rotational ehects are treated as 
small axisymmetric and stationary perturbations. The perturbed internal and external 
spacetimes are matched with a uniformly rotating thin shell. We assume that the 
angular velocity of the thin shell, hi, is much smaller than the Keplerian frequency 
of the nonrotating gravastar, The solutions within an accuracy up to the second 
order of ff/fffc are obtained. The thin shell matter is assumed to be described by a 
perfect fluid and to satisfy the dominant energy condition in the zero-rotation limit. In 
this study, we assume that the equation of state of the thin shell matter is given by 
dp/da = {dpo/dR)/{dao/dR). In other words, the equation of state for perturbations 
is the same as that of the unperturbed solution. The spherically symmetric component 
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of the energy density perturbations, 5 (To, is assumed to vanish independently of the 
rotation rate, i.e., Sco = 0 is always assumed. 

Using the nonrotating thin shell gravastar stable against radial perturbations 
derived by Visser and Wiltshire [3] as unperturbed solutions, under the assumptions 
mentioned before, we numerically evaluate the rotational corrections to the structure 
of the gravastar up to the second order of For nonrotating gravastar models, 

we choose hve sequences of equilibrium solutions characterized by hxed values of 
M/L = 1/7, l/\/40, 1/6, 1/5, and l/\/^- By changing values of R/M, for a given value 
of M/L, we may obtain one sequence of the non-rotating gravastar. The maximum and 
minimum values of R/M are determined for the thin shell to satisfy the dominant energy 
condition. For each sequence of the nonrotating gravastar models with M < 1/6, there 
is a singular non-rotating gravastar model in the sense that the master equations for the 
quadrupole perturbations have no solution. In other words, the gravastar cannot rotate 
perturbatively for such singular non-rotating gravastar models. These singular solutions 
have a value of R/M in the range of 2.9 < R/M < 3. To avoid the occurrence of such 
singular solutions, in this study, we focus on the nonrotating gravastar models with 
R < 2.8M and examine rotational effects on the structure of these gravastars. Note 
that some authors have shown that horizon-less compact objects with stable photon 
circular orbits are nonlinearly unstable [muH], which suggests that the gravstars with 
R < 2.8M are nonlinearly unstable. However, we think that because of the nonlinearity, 
it is required further investigations to say that this nonlinear instability is general, 
although it indeed operates for some situations. Interesting hndings in this study are 
the following: (i) Some rotating solutions have the dimensionless quadrupole moment, 
QM/J'^, the same as that of the Kerr black hole, i.e., QM/J^ = 1. This means that if 
those gravastars are incidentally observed, we cannot distinguish gravastars from black 
holes by using their values of the mass, M, spin angular momentum, J, and mass 
quadrupole moment, Q. (ii) Most gravastar models obtained in this study satisfy the 
relation of QM/J"^ < 1. This is absolutely different from that of the standard rotating 
compact stars like neutron stars, for which QM/ > 1 is satished. (iii) Some rotating 
gravastar models are prolate. This is again an opposite behavior to that in the standard 
rotating compact stars. The reasons for occurrence of (ii) and (iii) are not yet clear. 

The present slowly rotating gravastars become slowly rotating thin shells with 
isotropic pressure in vacuum asymptotically flat spacetime in the limit of the zero 
cosmological constant (in the limit of L —)■ cx)). Taking the limit of L —)■ ex, thus, 
we obtain the slowly rotating thin shell models and investigate their properties. For the 
thin shell, the dominant energy condition is satished if R/M > 25/12 2.08. Thus, we 

may take the Newton limit of the slowly rotating thin shell in vacuum asymptotically hat 
spacetime. We hnd that many rotating thin shells calculated in this study are prolate 
and that dimensionless quadrupole moment for most of them satisfy QM/J^ < 1 and 
become negative when the radius of the shell is larger than 2.9M. This means that under 
our assumptions, slowly rotating thin shells are prolate even in the Newton limit. Thus, 
we do not think that the appearance of the prolate-shape slowly rotating gravastar is 
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attributed to the general relativistic effects. 

We also examine properties of the ergoregion for the present thin shell gravastar 
models since its existence results in secular instability and their rotational energy can be 
extracted through the ergoregion instability. To examine the ergoreion of the thin shell 
gravastar, we basically follow the prescription devised by Schutz and Comins [26l |25] , 
in which the e^-order rotational effects are partly included in the spacetime for the 
sake of tractability. However, we also consider the case where the e^-order rotational 
effects are fully included in the spacetime. We hnd the following properties: (i) The 
ergoregion forms only outside the thin shell, (ii) Percentage differences between the 
minimum values of eJ/M^ at which the ergoregion first forms obtained with the Schutz 
and Comins’s treatment and by the spacetime including all the e^-order rotational effects 
are ~ 20%. (hi) This relative difference becomes smaller for more compact gravstar 
models as expected by Schutz and Comins |26l |25] . (iv) The ergoregion determined with 
the Schutz and Comins’s treatment is smaller than that determined by the spacetime 
including all the e^-order rotational effects. 
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Appendix A. Extrinsic curvature 


The extrinsic curvature is approximated by + e + O(e^) for 

slow rotation approximation. The explicit forms of the nonzero components of 
are given as follows: 
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Appendix B. The e^-order perturbations of the energy density and pressure 

The quadrupole components of the energy density and pressure perturbations are, 
respectively, given by 
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The spherically symmetric components of the energy density and pressure perturbations 
are, respectively, given by 
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The explicit form of B is given by 


where 


B = 


U 




^/J+ 3M-R\ J‘^{M + R) f 1 


R 


8Mi?3 


vr 


[M./J+{2J - R^Vtkf 2J^{M-R){QM^ 


R-M\ 

'rTFJ 

R^) 


2AR\f+f V 


/- 




(2J - R^Vtk){2J{3M -R)- R^{QM^ - R^)nk) 

7r 

R^VLk{R^ - MR - 3M^){2J + R^{R - 3M)fifc) 

7F 

2i? + 3M , L^/J^G LMW \ 

~WFf) 


AR^ 


+ 


^ / 

da \ 

M + R f {2J - R^Vtkf VLk{2J + R^{R + 3M)VLk) 


U 
R 

4Rf~ 


2F 


R^^/T 


LW 


J^{5M - AR) _ 

ARm^/J+^ ^ 3R^{f+YF 

J^(-8M^ + 9MR-3R^ 


X 


1 fR-3M 


M 

1 \ f 


M^{2J - R^^lkf 

+ MR^Vtk{2J + R^{R- 3M)Vtk) 


VM 


!^iv 

da 


-/rJ U^/F^/F ^ 

3M + 2R Ly/f^G LMW \ RLy/J^W 


^Rvr 


2F 


_^f 3R\R + 2M) ^^^^^ 


AM V 


2M 


rVTfJ 

QM^ + SM'^R 


F 

3MR^ - 


QliX) 


3R^ 


Rf+ 


-1 


(C.l) 


A = 


M 


1 , 


Y = 


R 


R 




2M’ 

R^Q:k - 2 J 


(C.2) 

(C.3) 


1 + 


MLf-G 


X 


R^n 


3f+ V Rf^R 
k-2J 2J{2R-3M) + R^{R-3M)VLk 


VT 


J2 


M/p 

1 
X 


1 / = 


77^ 

Lf-G 

FM 


R + 2M 

^ R^G 
1 + 


R^/r 

Lf-G{M + R) 


LF 


VT 


1 + 


MLf-G 

Rf-^F 


-1 


(C.4) 


2M^M+ 2R) + 3R^{R-3M) 3R^f 


2 f+ 


R^f^ 


+ 


2M 


logF 


















































Slowly rotating thin shell gravastars 


35 


2M2 - ?,MR - 3i?2 3i?(i?2 _ M^ogY | 


X 


W = 


MR 

1 RG 

r^Tf~ 


-3L‘^R 




1 + 


2M2 

MLf-G 

Rf^F 


-1 


2R'^ + 3L^f Arctanh(i?/L). 


(C.5) 

(C.6) 


References 

[1] Markov M A 1982 Limiting density of matter as a universal low of nature Sov. Phys. JETP Lett. 

36 265 

[2] Mazur P O and Mottola E 2001 Gravitational condensate stars Preprint gr-qc/0109035 

Mazur P O and Mottola E 2004 Gravitational vacuum condensate stars Proc. Natl. Acad. Sci. Ill 
9545 

[3] Visser M and Wiltshire D L 2004 Stable gravastars - an alternative to black holes? Class. Quantum 

Gray. 21 1135 

[4] Garter B M N 2005 Stable gravastars with generalized exteriors Class. Quantum Gray. 22 4551 

[5] Gattoen C, Eaber T and Visser M 2005 Gravastars must have anisotropic pressures Class. Quantum 

Gray. 22 4189 

[6] Ghirenti G B M H and Rezzolla L 2007 How to tell a gravastar from a black hole Class. Quantum 

Gray. 24 4191 

[7] Gardoso V, Pani P, Gadoni M and Gavaglia M 2008 Ergoregion instability of ultracompact 

astrophysical objects Phys. Rev. D 77 124044 

[8] Pani P, Berti E, Gardoso V, Ghen Y and Norte R 2009 Gravitational wave signatures of the absence 

of an event horizon: Nonradial oscillations of a thin-shell gravastar Phys. Rey. D 80 124047 

[9] Ryan E D 1995 Gravitational waves from the inspiral of a compact object into a massive, 

axisymmetric body with arbitrary multipole moments Phys. Rev. D 52 5707 

[10] Ghirenti G B M H and Rezzolla L 2008 Ergoregion instability in rotating gravastars Phys. Rev. D 

78 084011 

[11] Uchikata N and Yoshida S 2014 Slowly rotating regular black holes with a charged thin shell Phys. 

Rev. D 90 064042 

[12] Hartle J B 1967 Slowly Rotating Relativistic Stars. 1. Equations of Structure Astrophys. J. 150 

1005 

Hartle J B and Thorne K S 1968 Slowly rotating relativistic stars. 11. Models for neutron stars 
and supermassive stars Astrophys. J. 153 807 

[13] Thorne K S 1971 General relativity and cosmology edited by Sachs R K (Academic Press, New 

York) 

[14] Ghandrasekhar S and Miller J G 1974 On slowly rotating homogeneous masses in general relativity 

Mon. Not. R. Astron. Soc. 167 63 

[15] Israel W 1966 Singular hypersurfeces and thin shells in general relativity Nuovo Cimento 44B 1 

[16] Barrabes C and Israel W 1991 Thin shells in general relativity and cosmology: The lightlike limit 

Phys. Rev. D 43 1129 

[17] Keir J 2014 Slowly Decaying Waves on Spherically Symmetric Spacetimes and an Instability of 

Ultracompact Neutron Stars arXiv:1404.7036 

[18] Gardoso V, Crispino L C B, Macedo CEB, Okawa H and Pani P 2014 Light rings as observational 

evidence for event horizons: Long-lived modes, ergoregions and nonlinear instabilities of 
ultracompact objects Phys. Rev. D 90 044069 

[19] Brill D R and Cohen J M 1966 Rotating Masses and Their Effect on Inertial Erames Phys. Rev. 

143 1011 

Cohen J M 1967 Note on the Kerr Metric and Rotating Masses J. Math. Phys. 8 1477 







Slowly rotating thin shell gravastars 


36 


[20] de la Cruz V and Israel W 1968 Spinning Shell as a Source of the Kerr Metric Phys. Rev. 170 

1187 

[21] Pfister H and Braun K H 1985 Induction of correct centrifugal force in a rotating mass shell Class. 

Quantum Gray. 2 909 

[22] Pfister H and Braun K H 1986 A mass shell with flat interior cannot rotate rigidly Class. Quantum 

Gray. 3 335 

[23] Abramowicz M A and Prasanna A R 1990 Centrifugal-force reversal near a Schwarzschild black 

hole Mon. Not. R. Astron. Soc. 245 720 

[24] Friedman J L 1978 Ergosphere instability Commun. Math. Phys. 63 243 

[25] Comins N and Schutz B F 1978 On the ergoregion instability Proc. R. Soc. Land. A 364 211 

[26] Schutz B F and Comins N 1978 On the existence of ergoregions in rotating stars Mon. Not. R. 

Astron. Soc. 182 69 

[27] Yoshida S and Eriguchi E 1996 Ergoregion instability revisited - a new and general method for 

numerical analysis of stability Mon. Not. R. Astron. Soc. 282 580 

[28] Bardeen J M, Press W H and Teukolsky S A 1972 Rotating black holes: Locally nonrotating 

frames, energy extraction, and scalar synchrotron radiation Astrophys. J. 178 347 


